Externally-driven transformations of vortex textures in flat submicrometer magnets 
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Two effects of oscillatory transformations of vortex textures in flat nanomagnets due to the ap- 
plication of an external field or a spin-polarized electric current are analytically described with 
relevance to soft- magnetic structures of submicrometer sizes (whose thickness is significantly bigger 
than the magnetostatic exchange length). These are changes of a domain wall (DW) structure in 
a long magnetic stripe (oscillations between a transverse DW, a vortex DW, and an antivortex 
DW) and periodic vortex-core reversals in a circular magnetic dot which are accompanied by os- 
cillatory displacements of the vortex from the dot center. In nanostructures of smaller thicknesses 
(comparable to the exchange length), where nonlocal magnetostatic effects are very strong because 
of fast spatial variation of the magnetization, similar phenomena have been widely studied previ- 
ously. Here, the dynamics is investigated within a local approach including magnetostatic field via 
boundary conditions on solutions to the Landau-Lifshitz-Gilbert equation only. Both the DWs in 
stripes and vortex states of the dot are treated as fragments of a cross-tie DW. Despite similarity of 
the cyclic transformations of the ordering to the dynamics of more strongly confined nanomagnets, 
details of motion (trajectories) of the vortices and antivortices (Bloch lines) of the textures under 
study are different, which is related to prohibition of rapid jumps of the polarization of Bloch lines. 
In addition to the magnetization rotation about the direction of magnetic field or current polariza- 
tion, the evolution of textures is shown to relate to oscillatory changes of the direction of a cross-tie 
DW with respect to any arbitrary axis in the magnet plane accompanied by oscillations of the DW 
width. 

PACS numbers: 75.75.Jn, 75.78.Fg, 85.70.Kh, 85.75.-d 

Keywords: domain wall, Landau-Lifshitz equation, spin-transfer nano-oscillator, magnetic stripe, magnetic 
dot 
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I. INTRODUCTION 

Application of external (longitudinal) magnetic field 
or voltage to the DW-containing ferromagnetic nanowire 
drives translational motion of the DW. When the inten- 
sity of the field or electric current exceeds critical value 
of the Walker breakdown, the DW translation is accom- 
panied by cyclic transformations of the DW structure 1 -. 
In soft-magnetic nanostripes, strong magnetostatic field 
stabilizes textures (DWs) whose topological charges are 
pinned to the platelet boundary and whose topology gets 
transitions with changing thickness and width of the 
nanostripe 2 4 . There, the field-induced DW transforma- 
tions are transitions between the so-called transverse and 
vortex DWs (antivortex DWs) 5 8 . 

In circular magnetic dots, magnetostatic field can form 
a centered vortex of magnetization 9,10 . In the nanocon- 
tact with a spin valve structure, the dot conducts in- 
plane-polarized electric current in the out-of-plane direc- 
tion, which leads to the transition of vortex (with increase 
of the current intensity) to the state of precession about 
the center of the dotpii - — . At low current intensity, the 
vortex sustains displacement moving out the dot 15 . 

Electrically-induced dynamical transformations of pla- 
nar spin-textures play the key role in magnetic nano- 
oscillators for applications to microwave generation and 
sensing^"— . In terms of design of DW-based storage 
and logic devices^"—, the DW-structure transformation 
due to longitudinal magnetic field or electric current is 
an undesired effect which limits the efficiency range of 



velocities of the DW propagation 22 . 

Since the ordering in soft-magnetic nanostructures is 
governed by the exchange as well as by long-range dipo- 
lar interactions (permalloys are the most popular ferro- 
magnetic materials for the above applications), descrip- 
tion of dynamical phenomena is nonlocal and remains 
a challenge. At present, micromagnetic simulations are 
the main source of knowledge on details of the texture 
evolution in spatially confined magnets. Analytical stud- 
ies of the dynamics of DWs in nanostripes j 23 i 24 , and of 
vortex-states in circular dots^, use modifications of a 
Thiele method that treats the magnetic vortices as mov- 
ing rigid objects 2 ^. The spin-structure transformations 
are related to the motion of a number of vortices and an- 
tivortices, thus, the dynamical (Thiele) equations of the 
(anti)vortex-core positions constitute complex systems. 
Moreover, they contain semi-empirical (gyrotropic and 
viscosity) coefficients to be estimated. 

With relevance to thick-enough planar nanomagnets, 
in the present paper, I analytically study the field- and 
current-induced transformations of DWs in ferromag- 
netic stripes solving the Landau-Lifshitz-Gilbert (LLG) 
equation instead of using the Thiele equations. Within 
a relatively-simple local approach, the magnetostatic ef- 
fects are included via boundary conditions only. Also, I 
apply the present method to the vortex-state transforma- 
tions in a circular magnetic dot. I treat the vortex- state 
in such a quasi-2D system as well as the DW in the mag- 
netic stripe as fragments of an infinite spin texture called 
the cross-tie DW 27,28 , (recently this approach has been 
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applied to study DW interactions in Refi2S). 

In contrast to the most frequently investigated na- 
noelements of ultimately-small thicknesses, in (thicker) 
systems under considerations, the exchange interactions 
dominate over the magnetostatic effects while the later 
ones are included as a perturbation 3 -^. In very flat sys- 
tems, the dipolar interactions dominate while the ex- 
change is thought of as a perturbation^* 3 - 2 . Similar ap- 
proach is used by popular codes for micromagnetic sim- 
ulations that apply the fast Fourier transform to calcu- 
late the dipolar (magnetostatic) field, which limits their 
validity 2 ^ 33 . The crossover between both the interaction 
regimes is relevant to thickness of the magnetic element 
close to the critical length l c of a "macrospin approxi- 
mation" defined in Ref, 34 and estimated to be from the 
range Al ex -=- 8Z e£C , where l ex denotes magnetostatic ex- 
change length (l ex « 5nm for permalloys). In elements 
whose spatial sizes are smaller than l c the density of ex- 
change energy is too high to allow an inhomogeneous or- 
dering unless the anisotropy is nonuniform 3 -. Decrease of 
one of the spatial sizes below l c in presence of strong local 
anisotropy can result in discontinuity of spatial variation 
(differentials) of the magnetization. 

Solving the LLG equation for the thick-film regime, 
differently oriented static cross-tie DWs are found to be 
analogs of the transverse and vortex DWs of very thin 
stripes. Upon externally-induced driving, the DW tex- 
tures oscillatory transform between these two basic con- 
figurations, however, details of the transitions between 
transverse and vortex states are found to differ from 
those in the thin stripes in terms of trajectories of the 
DW-texture defects (vortices and antivortices). Unlike 
most of available descriptions of the vortex dynamics in 
nanodots which deal with defects in the in-plane-ordered 
(curling) state, here, evolution of a vortex texture dis- 
torted from the plane at the dot rim is analyzed. The 
trajectories of the current-induced vortex translation are 
examined as well as trajectories of antivortices which ap- 
pear due to enforced reversals of the defect core. 

In Sec. II, stationary DW solutions to the LLG equa- 
tion in the 2D stripe geometry are found. Externally- 
driven dynamics of such structures is analyzed in Sec. 
III. Sec. IV is devoted to study of the structure and dy- 
namics of vortex states in circular dots. Conclusions are 
collected in Sec. V. 



Here, M = |m|, J denotes the exchange constant, j3\ 
determines strength of the easy-axis anisotropy, H = 
(H x , 0, 0) represents the external (longitudinal) magnetic 
field, thus, 7 denotes the gyromagnetic factor. Al- 
though, in permalloy-like soft-magnetic materials, bulk 
anisotropy is negligible, for generality of considerations, 
I admit the presence of the easy-axis anisotropy (e.g. due 
to wire deposition on substrate^), however, assuming it 
to be weak compared to the exchange, f3\ <C J/w 2 , where 
w denotes the stripe width. In soft-magnetic stripes, the 
shape anisotropy due to dipolar interactions (vanishing of 
magnetostatic charges at the magnet surfaces) aligns the 
magnetization of homogeneously ordered domains onto 
the direction of long axis of the stripe i = (l,0,0)r»2Z. 
Therefore, I study solutions that satisfy the condition 
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D m = ±(M,0,0). 



Since (pQ) must be solved with the constraint on the 
magnetization length, it is comfortable to consider equiv- 
alent to (pQ) equations of the unconstrained dynamics. 



Introducing m± 



± im z , one represents the mag- 



netization components with a pair of complex functions 
g(x, z, t), f(x, z, t) (secondary dynamical variables); 



2M 



f*/g + g*/f 



n J g-g f , , 

m * = M 777 — , — r77> ( 2 ) 
f*/9 + 9*/f 



thus ensuring that |m| = M. Insertion of ([2]) into (pQ) 
leads, following the Hirota method for solving nonlinear 
differential equation ^ 38 ! 39 , to the trilinear equations of 
motion 

- fiD t f* -g = f [« A + J (Pi + A 2 )] f ■ 9 + Jg*(D 2 x 

+D 2 z )g-g-( 1 H x + p 1 )\f\ 2 g, 
-g*\D t f* -g = g* [aD t - J(D 2 X + D 2 Z )] f* ■ g - Jf(D 2 x 

+D 2 z )r-r + {- 1 H x +p 1 )\g\ 2 r, (3) 

where D tl D x , D z denote Hirota operators of differenti- 
ation 

D^b{x 1 z, t) • c(x, z, t) = 



(d/dx — d/dx ) n b(x,z,t)c(x ,z ,t )\ x= 



-X ^Z — Z ,t — t 



For H = 0, stationary single- DW solutions to ([3]) are of 
the form 



/ = 1, g = ue 



(4) 



II. DOMAIN- WALL STATES IN 
FERROMAGNETIC STRIPE 

Let us consider stationary DW solutions to the LLG 
equation in 2D 



dm 

~dt 



J_ 
M' 



-m x 



M 



d 2 m d 2 m 



dx 2 dz 2 



(m • z)m x i 



7m x H 
a dm 



(1) 



where 



7 2 , 2 Pi 

k +q =- 



(5) 



and Re/c 7^ 0. I denote k = k + ik , q = q + iq , where 
k ( ), q ( ) take real values, hence, (|5j) is equivalent to 



7 '2 , '2 7 "2 "2 Pi 

k z + q z - k z - q = —. 

k k + q q =0. 



(6) 
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Assuming one of the DW edges to be centered at x = 0, 
(then u = e 1( ^), the relevant magnetization profile [the 
single-DW solution to (pQ)] is written explicitly with 



a) 



m+(x,z) = Me iW+fc x+q z) s,ech[k x ■ 
m x (x, z) = — Mtanh[fc x + q z]. 



q z], 



(7a) 
(7b) 



Let q ^ since, in the opposite case, the DW states 
are similar to DWs in ID ferromagnet and cannot exist 
in absence of the bulk anisotropy while I consider soft- 
magnetic systems with very weak anisotropy admitting 
the case /3\ = 0, 39,40 . Defining = arctan(g jk ), via 
©, one finds k" = -q tan((9) and k' 2 -q 2 = /3i/{J[l + 
tan 2 (#)]}. Also, I assume the magnetization orderings 
on both the stripe edges to be similar, thus, the phase 
factor of the RHS of (|7a|) changes by nir along the DW 
line k x + q z = between it ends, where n = 1,2, — 
It leads to the condition k (—wq jk ) + q w = nir and, 
finally, to q = nn / {w[l + tan 2 (#)]}. 

Additional boundary condition is related to minimiza- 
tion of the surface (magnetostatic) energy and it discrim- 
inates between different values of </>, n, and 0. I evaluate 
the energy of the DW with dependence on these param- 
eters using the Hamiltonian H = Ho + Hz, where 



-( 
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dm 




+ A[ 


2M I 


dx 




dz 
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2M L 



M 



(m-if 



-7H • m, 



(8) 



(Hz denotes its Zeeman part). Total energy of the DW 
E = Eo + Ez + Eb is the sum of the bulk energy + 
E z , defined by £7 (z) = X!^ Jo H (z)^dx, and of the 
boundary energy Eb which is of the magnetostatic origin. 

Since evaluating the magnetostatic energy of the stripe 
within phenomenological approach is a complex math- 
ematical problem whose systematic solution is beyond 
the scope of this paper, I search for energy of the stripe 
boundary using a heuristic argumentation. The formula 
for Eb is determined referring to a theorem by Carbou 
who proved that the magnetostatic energy of any ferro- 
magnetic element of finite thickness r; ~ 1/A 2 J s (m • 

n) 2 ds tends to 1/A 2 J dS (m • n) 2 d/ with r O, 4 ^ 2 -. 
Here, S denotes the surface of the bulk ferromagnet and 
dS denotes the boundary of the base of its solid, n is nor- 
mal to the magnet surface, n denotes the unitary vector 
outward to the line of the base boundary. The coeffi- 
cient A2 scales with r and with a diameter w following 
A2 ~ A 2 /[t| log(r/^)|], (with relevance to a stripe, w rep- 
resents its width 31 ** 3 - 2 -). One has to notice that the Carbou 
theorem is not strictly applicable to systems with open 
boundaries, (infinite stripes), however, it shows some ten- 
dency in ordering at the stripe edges. In particular, it 
indicates that the magnetostatic interactions in fiat mag- 
nets induce more than one hard directions of magnetiza- 
tion parallel to the main plane whereas strength of the 
hard-axis anisotropy increases with size of the platelet 
along this axis. I propose to effectively describe the mag- 
netostatic energy in the form of an integral over the stripe 
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FIG. 1. DW configurations; a) a transverse DW, b) a vor- 
tex DW. In the upper draws, arrows indicate magnetization 
alignment. 



edge 



E B 



2 \ ^ 2 

■™>x) + T- m * 
A 2 



dx (9) 



while e.g. the approach of Ref^ 2 - corresponds to 1/Ai = 
0. By analogy to finite platelets, this coefficient is ex- 
pected to scale with the stripe width following Ai oc 
X 2 /w. More detailed estimation of Eb is performed in 
Appendix A. 

I mention that if x direction was not one of the local 
hard axes at the boundaries of a wire (stripe) of finite 
length, the DW would be unstable due to a magnetostatic 
field at the wire (stripe) ends 43 . This local anisotropy 
prevents spontaneous DW motion towards one of the wire 
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ends. 

Inserting 
rives at 



a|)- ([7b|) into the Hamiltonian (|8j), one ar- 



E o (0,n) = 2JMw 



tan 2 



(10) 



Evaluating Eb, I divide it into two parts (Eb 

Eb2)] 



Ebi 



E B1 {0,n) = — / sech <^ — f dx, 

Ai [w[l + tan 2 ((9)] J 



E B2 (^0,n) 



2M 

T9. 



2 />oo ( 

/ sin 2 , 



n7r tan(#)x 
^[1 + tan 2 (<9) 



xsech 



;[1 + tan 2 



dx. 



(ii) 



The minimization of E B {(j),0,n) leads to <p = 0,7r inde- 
pendently of other parameters of the DW. These values 
of (j> correspond to the presence of halfvortices (halfan- 
tivortices) at the stripe edges (as shown in Fig. 1). In 
the regime of narrow stripes w/r ~ 1, A2 <C Ai, (I/A2 
is small, thus, Eq dominates over Eb), via minimiza- 
tion of Eo(0,n) the smallest possible value of n is prefer- 
able while the minimization of both Eq{0, n), £^#(0, n) 
points out = 0, 7r to be preferred. Increase of the stripe 
width w with fixed thickness r results in transition to the 
regime A2 > Ai. Furthermore, since j3\ <C J/w 2 , from 
([TDD , Eo(6,n) « Eo(n) is independent of w, while, from 
(fTTjh £"5 oc w/Ai. Hence, for big-enough w, Eb becomes 
comparable to Eq. By the minimization of EBi(0,ri), 
the biggest possible value of tan 2 (#) and the smallest 
value of n are preferable, whereas, for ^ 0, 7r, the min- 
imization of i?B2(0,#,n) indicates big values of n to be 
preferable. Hence, the transition from the DW state of 
n = 1, 6 = 0,7r to a state of n = 2, 7^ takes place 
with increase of w. For n = 2, the condition tan 2 (#) > 1 
has to be satisfied in order that E B (0, 9, 2) < Eb(0, 0, 1). 

The state of n = 1 and # = , 7r correspond s to 
/c" = <?' = 0, \q\ = tt/w, \k'\ = ^tt 2 /w 2 + f3i/J and 
it is called a transverse DW. With relevance to the state 
of n = 2, I take /?i = for simplicity, and \0\ to be close 
to its infimum \0\ = 7r/4. The resulting magnetization 
structure corresponds to \q \ = \q \ = \k \ = \k \ = ir/w 
and one calls it a vortex DW. The polarity of vortex 
(transverse) DW, (the magnetization orientation in the 
center of vortex/halfvortex, parallel or antiparallel one 
to y axis), is determined by value of <p while q /k = ±1 
determines its chirality, (the direction of magnetization 
rotation in the stripe plane in the vortex or antivortex 
cores, clockwise or anticlockwise one). 

Treatment of the ordering of soft-magnetic stripe 
within the present local approximation of micromagnet- 
ics provides a unified description of transverse and vor- 
tex DWs which is consistent with previously known exact 
DW solution for the so-called exchange dominated regime 
of ultra-thin nanostripes (n = l) 31 i 32 , and for 2D stripes 
with exchange interactions of Ref. 44 , (n = 2). 



III. DYNAMICAL TRANSFORMATIONS OF 
DOMAIN WALL 

The application of a longitudinal magnetic field (H ^ 
0) enforces the DW translation along the wire (stripe) 
whose direction (parallel or antiparallel to x axis) is de- 
termined by the decrease of Zeeman energy with time. 
Two regimes of the field intensity, separated by a crit- 
ical (Walker breakdown) value Hw, have to be distin- 
guished. For weak field \H X \ < ELw, the dynamics is 
restricted to the stationary DW translation, whereas, for 
\H X \ > Hw, the magnetization rotation about the long 
axis of the stripe stripe is allowed. The phenomenon 
of Walker breakdown has been predicted by solving the 
LLG equation in the framework of ID model with two- 
axis anisotropy 1 . Although, numerical estimations of 
Hw for nanostripes show invalidity of the ID model to 
the systems under study, the field- and current-induced 
Walker breakdowns are observed in them 4 ^. 

In order to include the breakdown phenomenon into 
the description of flat magnets, I discriminate between 
"weakly dissipative" and "purely relaxational" dynami- 
cal regimes. By analogy to the construction of a time- 
dependent Ginzburg-Landau equation in the theory of 
dynamic critical phenomena, in the "weakly dissipative" 
regime, the evolution is governed by Eq. (pQ) while the 
"purely relaxational" dynamics corresponds to a modi- 
fied [by changing the LHS of (pQ) into zero] equation of 
motion. I claim the conservative kinetic term (the LHS) 
of (pQ) to be irrelevant to the case of a weak-field (-current) 
since such a field induces a low-energy spin excitations 
that have to be overdamped in an anisotropic medium 4 ^. 
In the stripe model of Sec. II, the anisotropy is local (it is 
introduced via boundary conditions) and the mechanism 
of breakdown accompanied by disappearance of the ki- 
netic term in (pQ) is similar to overdamping low-energy 
excitations (the central peak in excitation spectra) of 
isotropic systems doped with anisotropy centers 47,48 . 

Searching for the field-driven DW evolution, I modify 
the ansatz (j4]) into 



/ = i, 



g = ue kx+qz - lt . 



(12) 



With relevance to the "weakly dissipative" dynamical 
regime of \H X \ > Hw, via (|3]), for k and q of the pre- 
vious section, I find 



I = 



-lH x 



a 



(13) 



The imaginary part of I is equal to the frequency of the 
magnetization rotation about x axis. This rotation will 
be shown to imply additional oscillations of the orienta- 
tion of DW in the stripe plane due to effect of the mag- 
netostatic field at the stripe boundary. Below the break- 
down, for \H X \ < Hw, the DW solution to the "purely 
relaxational" secondary evolution equation [the LHSs of 
(j3j) are changed into zero] corresponds to I = jH x /a. In 
both the cases of "weakly dissipative" and "purely relax- 
ational" dynamics, motion of the cross-tie DWs can be 
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thought of as a translation of vortices and antivortices of 
the DW texture. For \H X \ < Hw the velocity of this mo- 
tion is constant and oriented along the stripe, whereas, 
for \H X \ > Hw, its direction is not conserved. In the last 
case, since directions of the vortex (antivortex) transla- 
tion deviate from x axis, value of the energy Eb evolves, 
which results in an evolution of the DW parameters 0, n. 
In order to determine 0(t), n(t), I minimize the energy 
£ [0(t) , n(t)} + E B [jH x t/(l + a 2 ) , 0(t) , n(t)) with respect 
to these functions. Since I assume 6(t), n(t) to be inde- 
pendent of spatial variables, the DW stays in a cross-tie 
structure during whole the evolution time, which is due 
to domination of the exchange interactions over the dipo- 
lar ones. Hence, the density of Bloch lines (the distance 
between vortices) and the DW width are uniform along 
the wall^, however, unlike in theories of the Bloch wall 
and Bloch line motion in platelets (which are valid to 
small propagation distances)^, they are not conserved 
as far as \H X \ > Hw- 
The constraint 

n(£)/{l + tan 2 [#(£)]} = l (14) 

corresponds to transitions between states (n = 1,0 = 
or = 7r), (n = 2, \0\ = 7r/4), . . ., transverse DWs and 
(multi-) vortex DWs. The motivation for writing ([T4]) is 
following. Along any straight z = zo G (0, w), the magne- 
tization component m x (x,zo,i) sustains a shift without 
changing its profile, {m x (x, zo, t) = m x [x — xo(t), zo, 0]}. 
It is because the applied longitudinal field H drives 
only rotation of the remaining magnetization compo- 
nents about x axis while any straight z = zq does not 
intersect the stripe edges where a local anisotropy is 
present. This condition is fulfilled when the DW pa- 
rameter k = ±q is conserved and q is transformed into 
q (t) = =F& (t) = k tan[0(t)], which is equivalent to ([14]) . 
Inserting the unperturbed magnetization components 

m+(x, z, t) = Me [ (- l '' t+k '' x+q " z) sech[k x + qz - I't], 
m x (x, z, t) = — Mtanh[/c x + q z — I t], (15) 

(with l' = Re/, l" = ImZ) into the sum Eo + Eb, for 
f3i « and \l 1 1 | = a <C 1, and applying (fl4|) . one finds 

Eb(a) = 2JM7r(l + a 2 ), 

„ , 2M 2 Z 400 2 /7TX\ , 
E B (a,t) = — / sech — dx 

+ / S in 2 7 ^V1 + a 2 ) 

A2 J_oo 1 ^ J 

xsech 2 f— )dx, (16) 

where, a = tan(^). The variable a is a genuine variational 
parameter for the present extremum problem, (a = cor- 
responds to the transverse DWs, a = ±1 to the vortex 
DWs). For some arbitrary chosen values of J, Ai, A 2 , 
w, the energy function Eq(o) + Eb(cl, t) is plotted in Fig. 
2a. The corresponding density plot of the derivative of 
Eo(a) + E B (a : t) over a is presented in Fig. 2b where 




h ' 
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a 

FIG. 2. a) Sum of the exchange and boundary energies 
Eo + Eb with dependence on DW direction and time for ar- 
bitrary chosen values of the parameters JA\/(Mw) = 0, 026, 
JA2 / (Mw) = 0, 020. b) The corresponding density plot of the 
derivative of Eo(a) + ^s(a, t) over a. The contour in the plot 
center represents the solution to d[Eo(a) + EB(a,t)]/da = 0. 



the contour d[Eo(a) + E B (a,t)]/da = is shown. This 
contour indicates three extremal curves a(t) which relate 
to different trajectories of vortices and antivortices of the 
DW structure. The straight a(t) =0 does not minimize 
Eo(a) + EB(a,t), as seen from Fig. 2a, while the mag- 
netization rotation about a constant direction with the 
frequency \l | = j\H x \/(l + a 2 ) corresponds to a{t) > 
or a(t) < 0. The trajectories of vortices and antivortices 
in the DW texture are sketched in Fig. 3, where peri- 
odically created and annihilated transverse and vortex 
(antivortex) DWs are indicated with thin solid lines. 

It is seen from Fib)) that the dependent on A 1 contribu- 
tion to the energy of the DW [E B \ of (G2P/ does not play 
any role in the dynamical transformation of the texture 
since it is neither dependent on the variational parameter 
a nor on time. 



IV. VORTEX IN CIRCULAR DOT AND ITS 
TRANSFORMATIONS 

I study transformations of a magnetic vortex in a 
relatively-thick dot (of thickness r > 4Z e£C ) in a pillar 
nanocontact structure that contains a spin valve. The 
application of voltage to such a structure induces spin- 
polarized electric current through the dot plane and en- 
forces cyclic magnetization motion. Unlike in very thin 
dots, where vortices are of small radii compared to dot 
radii because of strong effect of the magnetostatic field^, 
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FIG. 3. Trajectories of vortices (dashed line) and antivortices 
(dash-dotted line) in the field-induced motion of a head-to- 
head DW. The arrows indicate the magnetization directions 
in relevant DW areas. 



in thicker exchange-ordered dots, the radius of vortex 
core is comparable to the dot radius and the out-of-plane 
component of the magnetization at the dot boundary can 
be nonzero52. 

Since skyrmions are static solutions to the isotropic 
LLG equation in 2D—, they seem to be natural can- 
didates to represent the vortices in dots. Indeed, in 
strongly flattened nanodots, ordering of vortex cores cor- 
responds to skyrmion profiles (see Ref*5^ for review of 
models of the centered vortex in the nanodot). Analytical 
approaches to the vortex- state dynamics of Refs j 25 i 55 are 
based on transformations between local two-skyrmion or 
skyrmion- ant iskyrmion textures embedded in an in-plane 
ordered curling background. They follow the observa- 
tion that vortex in such a state can be shifted from the 
dot center and deformed due to interaction with another 
(virtual) vortex without production of surface magneto- 
static charges 56 . However, in circular dots thick enough 
to allow using local approximation of the magnetostat- 
ics, the centered skyrmion is found not to be a stable 
state, (Appendix B), while any significant shift of the po- 
sition of static vortex from the dot center is not confirmed 
experimentally 57 . Moreover, when the dot is built into 
a spin- valve-containing pillar structure, a weak magnetic 
anisotropy can be present whose easy axis is directed in 
the dot plane. In the presence of such an anisotropy, 
skyrmions are not solutions to the LLG equation while 
cross-tie DWs are. Therefore, I consider a single Bloch 
line of the cross-tie DW to be the center of a vortex in the 
magnetic dot, which is consistent with an observed break- 
down of rotational symmetry of the dot orderin g 57 ! 58 . 

Studying magnetization configurations in the circular 
dots for H = 0, I use a static solution to (pQ) in the 
form (|Ta) - (|7bJ) with k" -q" q /k' -g"tan(0), k' 2 - 
q" 2 = Pi/Ul 1 + tan 2 (/9)]} « 0. Following the analysis of 
the previous section, the assumption of the orderings at 
the ends of any dot diagonal to be similar to each other 
leads to the relations q = nir / '{2R[1 + tan 2 (#)]} = ±k , 
where n = 1, 2, . . ., where R denotes the dot radius. The 



exchange energy of such a dot state takes the form 

r 27r r R c / 

E (a) = JMk 2 (1 + a 2 ) J J sech 2 |jk r 

x [sin((/?) + a cos ((/?)] |rdrd(/?, (17) 

where a = tan(0). Specific states characterized by <p — 
or <j) = 7r correspond to the presence of vortex or an- 
tivortex in the dot center. The boundary energy (of the 
magnetostatic origin) 



E B = 1/A 2 f (m-n') 2 dZ 

JdS 



(18) 



(where n denotes the unitary vector outward to the dot 
boundary and A 2 ~ A 2 /[r| log(r/2,R)|] scales with the 
dot thickness r) is evaluated with 

M 2 R f 2lT / ( i ~i 

Eb(cl, (j>) = — ^ — J y — tanh |/c R[sm(ip) + a cos ((/?)] j 

x sm(<p) + sech j/c R[sm(cp) + acos(y?)]| 
x sin |(/> + k R[cos(<p) — a sm(ip)] j 

M)V (19) 



X cos 



The plot of energy Eo(a) + Eb(cl,0) corresponds to the 
cross section of 3D plot in Fig. 4a at t = 0. Looking at 
this this cross section, one sees that absolute energy min- 
imum corresponds to|a|«l(|0|«7r/4) which represents 
a vortex DW (Fig. 5a). 

One includes an electric current perpendicular to the 
dot plane and spin-polarized in this plane via adding the 
spin-transfer torque 

— m x (m x p) 

M v ' 

to the RHS of (pQ). Here a is proportional to the cur- 
rent intensity and p denotes its spin polarization (|p| = 
1)25*55*52. Applying the transformation ©, for f)i « 
and p = (1,0,0), one arrives at the secondary equations 
of motion 

- fiD t f* -9 = f [aD t + J[Dl + D 2 Z )] f ■ 9 
+Jg*(D 2 x + D 2 z )g-g-ia\f\ 2 g, 

-g'iDtf -g = g* [+aD t - J(D 2 X + D 2 Z )] /* • g 

-Jf{D 2 x + D 2 z )r •/* -H«7| 2 r.(20) 

These equations can be obtained from (|3]) via change of 
jH x into icr, thus, insertion of the ansatz (fT2j) into ([20]) 
leads [by analogy to ([T3]) ] to 



I 



-1 + ia 



(21) 



In the weak current regime, below the Walker breakdown 
a < <Jw, the LHSs of ([20]) should be changed into zero 
and then 



I 



-1(7 



il . 



(22) 
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FIG. 4. a) Sum of the exchange and boundary energies 
Eo + Eb with dependence on DW direction and time for ar- 
bitrary chosen values of the parameters J, R, and A2 = A, 
[J A/ (MR) = 0,022]. b) The corresponding density plot of 
the derivative of Eo(a) + Eb(cl, t) over a. The contour in the 
plot center represents the solution to d[Eo(a)+EB(a>, t)]/da — 
0. c) The corresponding plot for A2 = 2A, and d), for 
A 2 4A. 



Since, in this regime, the real part of I is equal to zero, 
unlike magnetic field, the polarized current does not in- 
duce any translation of the DW while it is responsible 
for rotation of the magnetization in the wall area about 
x axis. 

Including the magnetostatic field at the dot boundary, 
I analyze vortex-state dynamics as an evolution between 
different cross-tie DWs, thus, following the method of 
previous section, I look for the time dependence of the 
DW parameters 0, n. I use the constraint ([T4]) which is 
motivated by considering the m x magnetization compo- 
nent along the straight z = 0. For a < aw, this compo- 
nent is conserved [m x (x, 0,t) = ra^x, 0,0)] upon appli- 
cation of a longitudinally-polarized current because this 
current drives the dynamics of the remaining magnetiza- 
tion components only, thus, it does not induce deviation 
of the magnetization from the local easy plane yz at the 
ends of the dot diagonal z = [at (x,z) = (±i?, 0)]. 
On lines z = zo ^ 0, the profile of m x is also conserved 
while its center sustains a shift {m x (x, zo,i) = m x [x — 
xo (t) , zo , 0] } because the current induces a deviation of 
m from an easy plane at (x, z) = (±\/i? 2 — Zq, zo). With 
the constraint (fl4]h transforming the texture parame- 
ter q into q (t) = =pfc (t) = k tan[0(t)], I substitute <\> 
in (fT9j) with I t and plot the energy Eo(a) + ^(a,t) 
for some arbitrary values of J, A2, R (Fig. 4a) as 
well as its derivative over a (Figs. 4b-4d). The con- 



z/2R z/2R 

, . . . \ _o_ 5 L^^^^J -0.5 f . . . . 

0.5 5 -0.5 0.5 -0.5 0.5 

x/2R x/2R x/2R 








0- 



FIG. 5. In picture a), density plots of magnetization compo- 
nents of the dot state \a\ = 1, the grey scale is linear in the 
range [-1,1]. Below, schemes of consecutive textures (of the 
head-to-head DW type) and trajectories of vortices (dashed 
lines) and antivortices (dash-dotted lines) during the transfor- 
mation induced by a spin-polarized current of the magnetic 
dot; b), d) A2 < A c , c), e) A2 > A c . The arrows in b) and c) 
indicate the magnetization alignment in relevant dot areas. 




FIG. 6. Consecutive textures (of the head-to- head DW type) 
in a strongly driven (a > aw) circular dot. The arrows indi- 
cate the magnetization directions in relevant dot areas. 



tour d[Eo(a) + EB(a,t)]/da = indicates the trajectory 
a(t) = tan[0(t)]. There exists a critical value A c such that 
the function a(t) is discountinuous for A2 > A c . Jumps 
of a(t) are connected to changes of its sign as follows 
from Fig. 4d. For A2 < A c , sign of a(t) is constant (see 
Figs. 4b, 4c). A number of consecutive textures created 
during a single period of the vortex-structure transfor- 
mation is visualized in Figs. 5b, 5d together with shape 
of the corresponding vortex and antivortex trajectories 
(Figs. 5c, 5e). 
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For a current of high intensity a > aw, the rela- 
tion ([2T]) leads to prediction of a texture motion along 
the polarization direction p = (1,0,0) with the velocity 
v = Rel/k = a/[k (1 + a 2 )]. In the confined geometry, 
however, this motion is suppressed when a decrease of en- 
ergy due to the spin-transfer-induced shift of the vortex 
core from the dot center equals an increase of the bound- 
ary energy Eb- A non-zero value of I corresponds to 
the magnetization rotation in the DW area (in the vortex 
core) with frequency I = a a /(I + a 2 ) about x axis. In 
order that such a rotation would not induce any instabil- 
ity of the structure (the boundary and exchange energies 
were constant), it should be accompanied by a preces- 
sion of the vortex (the DW texture) in plane of the dot 
shown in Fig. 6. The above determined frequency I of 
the vortex precession in a circular dot coincides with the 
relevant formula obtained in Ref. 60 by solving the Thiele 
equation where it has been verified with micromagnetic 
simulations. The radius of the vortex precession can be 
small compared to the dot diamete r 57 ! 58 . 



V. CONCLUSIONS 

Cyclic transformations between transverse and vortex 
DWs in a strong longitudinal magnetic field have been 
described with relevance to magnetic stripes with dom- 
inance of the exchange interactions over magnetostatic 
ones. Similar effect has been observed in experiments and 
simulations for very thin magnetic stripes (of the thick- 
ness of up to four magnetostatic exchange lengths, which 
corresponds to 20nm for permalloys), where magneto- 
static field is strong in bulk and dynamical distortions 
of the DWs take place. Such DW distortions cannot be 
included into the present approach. The predicted dy- 
namics can be verified with thicker stripes whose DW 
structures are expected to be more robust. 

A specific feature of the DW dynamics above the 
Walker breakdown are consecutive pinnings and depin- 
nings of a single halfvortex or an halfantivortices to the 
alternating stripe edges during the transformations of 
transverse DW into vortex DW and vice versa 2 ^. It is 
similar to the behavior shown in Fig. 3, however, since 
there is no any DW distortion in the present study, the 
half(anti)vortices are not strongly pinned to one of the 
stripe edges while they move away from the stripe on a 
small distance. Another difference with respect to DW 
dynamics in very thin nanostripes is lack of a simulated 
and explained with Thiele equation effect of dynamical 
reversal of the vortex-core polarity during its collision 
with the stripe edge« i 24 i 61 . It is because of different 
structures of the " magnetostatically-ordered" DWs of the 
cited works compared to the ones of "exchange-ordered" 
walls. In thin nanostripes, the field-excited texture can 
take the form of a transverse DW of total skyrmion num- 
ber G = 0, (G = ^2 a p a Qa with p a = ±1 denoting the po- 
larity and q a = ±1/2, ±1 denoting the chirality of ath de- 
fect in the DW texture 62 ), or a vortex DW of G = ±1. It 



means that the polarities of halfvortex and halfantivortex 
in the transverse wall could be the same while the polar- 
ities of halfantivortices at the edges of a vortex DW can 
be opposite. In the present considerations, only values 
of G — ±1 for the transverse DWs and G = ±2 for the 
vortex DWs are allowed, (the DWs are of a rigid cross-tie 
structure) . 

Two intensity regimes of the spin-polarized electric 
current are found to relate to periodic transformations 
of the vortex states in circular dots with frequencies of 
different orders of magnitude. Since the increase 

of the oscillation frequency I with a is much faster for 
currents weaker than a threshold value a < aw, (when 
I = a/a), than for stronger currents a > aw, [when 
I = aa/(l + a 2 )]. The regime a > aw corresponds to 
a vortex gyration in the dot plane, (this regime of the 
current intensities is used to be studied with relevance 
to the dc-induced generation of microwaves), while the 
regime a < aw corresponds to a vortex motion along 
an open trajectory. Most frequently, the vortex motion 
has been studied dealing with very thin dots with cur- 
rent above the threshold intensity. Schematic snapshots 
of the gyrating vortex of Fig. 6 can be compared to the 
ones observed in relatively thick (30nm) and of big diam- 
eter (3/im) dots during the magnetic- field-induced gyra- 
tion [since descriptions of both the processes are similar 
except change of the parameter (|2T|) into ([T3]) ] 58 . Small 
number of data on the vortex dynamics below the thresh- 
old is available. Similar (non-gyrotropic) evolution of the 
vortex state to the one of Fig. 5a has been observed in 
hexagonal dots (of 50nm thickness) upon application of a 
nanosecond electromagnetic pulse (such a pulse induces a 
spin-polarized current and a longitudinal magnetic field 
simultaneously) 63 . 
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Appendix A: Estimation of boundary energy 

The magnetostatic energy of a magnetic element con- 
tains contributions that relate to interactions of surface 
charges, volume charges, and interaction between surface 
and volume ones 

E MS = J J ^^dV(x)dnx') 

+ / / ¥^ d ~') 

+ / / ^^dS(x)dV(x'), (Al) 
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where p = — V • m, a = n m. Influence of magnetostatic 
interactions on the magnetization of the main body of 
any ferromagnetic platelet decreases with increase of its 
thickness, while it is not so in terms of magnetization 
of the platelet boundary. Following Ref. 64 , reducing one 
of the spatial dimensions with relevance to flat systems 
of thickness r and neglecting volume and base-surface 
terms, the above expression is transformed into the en- 
ergy of the boundary of 2D system (up to the multiplayer 
r) 

tE b =t 2 [ [ <r(x)<r(x')ln(|x-x'|/r)dZ(x)dZ(x') 
Jos J as 



I I 

JdS JS ha 



cr(x)p(x ) + p(x)cr(x ) 



xln(|x-x |/r)d/(x)d5(x 



(A2) 



Here Sbase denotes the surface of the platelet base. 

For a circular dot, p(x) is not defined at its boundary, 
thus, I neglect the second term and obtain 



tE 



3 =t 2 [ [ (n-m)(x)(n-m)(x) 

JdS JdS 

ln(|x-x / |/T)d/(x)d/(x) 

27rRr 2 \n(R/r) [ (n • m) 2 (x)d/(x). 
JdS 



(A3) 



For any DW in the stripe, p(x) ~ —dm x /dx ~ 
[M 2 - ml] /(MS) with a constant S close to DW 
width, {since —dm x /dx = M d[tanh(;r / 5 + cz/S)]/dx — 
[M 2 - M 2 tanh 2 (x/£ + cz/6)]/(M6), and \dm z /dz\ < 
\dm x /dx\}. I estimate the boundary energy with 

tE b =t 2 ( [ (n-m)(x)(n-m)(x , )ln(|x-x / |/r) 

JdS JdS 

(n-m)(x)— -(x ) 



xd2(x)d^(x ) -r 



/ / 

JdS Js ba 



dx 



drrir 



dx 

25r 2 \n(S/r) 



(x)(n-m)(x ) 



ln(|x-x |/r)dZ(x)d5(x ) 



ra 2 (x, 0, 0)dx 



/oo 
[M 2 — m 2 x (x, 0, 0)] dx, (A4) 
-oo 

where nw corresponds to an effective thickness of the 
surface layer of the stripe edge over which the magne- 
tization is independent of normal coordinate 2;, (k « 1 
and k oc 8). Since 5 ~ w, one arrives at Eb of ([9]). 

In order that boundary conditions to the LLG equation 
(within 2D approximation) were determined, one requires 
the relation l 2 x /w 2 <C r/wln(w/r) to be satisfied follow- 
ing Kurzke^, while, unlike it is assumed in his thin- film- 
limit approach, another condition r/w <C l 2 x /w 2 is not 
fulfilled in present considerations. In Ref. 42 , the last in- 
equality was used in order to show that the bulk-bulk and 
bulk-boundary terms of the magnetostatic energy vanish 
faster than the boundary-boundary term in the thin-film 



limit, thus, they are negligible for very thin structures, 
(also, then, a strong easy-plane anisotropy justifies using 
XY model). DWs of this "exchange" limit are described 
with exact formula in Refs i 31 i 32 . When r > l ex (the 
platelet is thick-enough), the bulk-boundary term can be 
non- negligible. Its importance grows with increase of lat- 
eral diameter of the platelet [stripe width, increase of the 
coefficient n of (|A4|) ] with constant r. In most of experi- 
ments (and in the present paper), sizes of the system lie 
far from the Kurzke regime, approaching to the region 
r/w^l 2 ex /w 2 . 



Appendix B: Instability of skyrmion in circular dot 

I determine the energy of a skyrmion centered with 
respect to a circular magnetic dot using a solution to 
(PQ) for P\ = 0, = (in absence of the anisotropy 
and external field), with the boundary conditions; m y — » 
— M, m z + im x — )► with \/ x 2 + z 2 — )► 00. Applying 
the Hirota method, via transforming magnetization to 
secondary dynamical variables 

gf 



m. 



1TO, 



2M 



12 , |„|2' m V 



I7F+W 



the skyrmion solution to the resulting equation ((3]) is 
found to be of the form 



/ = i, 



3 i77+ig-arctan(:r / z) 



(B2) 



where 7Z is a characteristic texture width, q = 1 
(skyrmion) or q = — 1 (antiskyrmion). The exchange and 
magnetostatic energies of the skyrmion and antiskyrmion 
are evaluated with flS) and ([T8]) 



E (TZ) = 



JM 



2tt r R 



= JMlT 



2/T Z 2 + 16r 2 /K 4 + 2r 4 /^ e 
(1 + r 2 /^ 2 ) 4 

r 2 /tz 2 + sr^/h 4 + 2R 6 /n 6 



{i + R 2 /n 2 y 



rdrd(p 
(B3) 



M 2 

E B ('7l,r ] ) = — 
A2 



2tt 



R 2 /1Z 2 cos 2 [(l - q)p - rj] 







(l + R 2 /K 2 ) 2 



Rdtp. 



[1^2)2 ^ cos 2 (77) g = l, 



AL 

A, 

a 2 (i+R 2 /n 2 ) 2 2 



q=-l. 



(B4) 



Via minimization of Eo(7Z)-\-Eb(7Z, rj), infinitely big tex- 
ture radius 7Z is preferred, which leads to the instability 
of the centered skyrmion (antiskyrmion) state. 

It should be emphasized that the above conclusion is 
valid as far as the local approach is suitable. Full eval- 
uation of the magnetostatic energy can indicate stabil- 
ity of a centered vortex embedded in a curling in-plane 
magnetization state, especially for dots of high radius-to- 
thickness ratic 54 . Such stable centered vortices are often 
called skyrmions because of skyrmion-like core magneti- 
zation, however, in a nomenclature widely used in differ- 
ent branches of physics, they are not skyrmions since they 
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are not defects in the ferromagnetically ordered back- ground. 
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